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MASS-DEPENDENT BARYON ACOUSTIC OSCILLATION SIGNAL AND HALO BIAS 
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ABSTRACT 

We characterize the baryon acoustic oscillations (BAO) feature in halo two-point statistics using 
iV-body simulations. We find that nonlinear damping of the BAO signal is less severe for halos in the 
mass range we investigate than for dark matter. The amount of damping depends weakly on the halo 
mass. The correlation functions show a mass-dependent drop of the halo clustering bias below roughly 
90 ft. _1 Mpc, which coincides with the scale of the BAO trough. The drop of bias is 4% for halos with 
mass M > 10 14 /i _1 Mq and reduces to roughly 2% for halos with mass M > 10 13 /i _1 M . In contrast, 
halo biases in simulations without BAO change more smoothly around 90 ft, Mpc. In Fourier space, 
the bias of M > 10 14 H~ 1 Mq halos decreases smoothly by 11% from wavenumber k = 0.012 ftMpc -1 
to 0.2 /iMpc -1 , whereas that of M > 10 13 /i -1 M Q halos decreases by less than 4% over the same 
range. By comparing the halo biases in pairs of otherwise identical simulations, one with and the 
other without BAO, we also observe a modulation of the halo bias. These results suggest that precise 
calibrations of the mass-dependent BAO signal and scale-dependent bias on large scales would be 
needed for interpreting precise measurements of the two-point statistics of clusters or massive galaxies 
in the future. 

Subject headings: Cosmology: theory — large-scale structure of Universe — dark matter 



1. INTRODUCTION 



Baryon Acoustic Oscillati ons ( BAO, iPeebles fc Yul 
Il970t iSunvaev fc Zeldovichl [19701) produce an im- 
print on the matter distribution by the perturba- 
tions in the photon-baryon fluid before recombina- 
tion. It pr ovides a standard ruler for probing dark 
energy Ce.g.. lEisenstein fc Hul 119981 iCoorav et all 120011: 
Blake fc Glazebrookl 120031: IHu fc Haimanl 120031: iLinderl 
20031: iSeo fc Eisensteinl 120031 ). and" is considered to 



be least affected by systematics ([Albrecht et al.l [20061 ). 
Thus far, the BAO signal has been detected in sam- 
ples f rom Sloan Digital Sky Su r vey (lE isenstei n et all 
120051: IPadmanabhan et al.l 120071: iPercival et all I2010D . 
2dF Galaxy R edshift Survey (iCqle et al.l I2005D. 6dF 
Galaxy Surv ey (iBeutler et all 120111) , WiggleZ Dark En- 
ergy Survey ([Blake et al.112011), and Baryon Oscillation 
Spectroscopic Survey (Anderso n et al.l I2012D . Future 
spectroscopic and imaging surveys, such as BigBOSS, 
Euclid, and Large Synoptic Survey Telescope, all aim 
to determine cosmic distances to sub- percent level pre- 
cision with the BA O technique (e.g.. [Zhan et al.l l2009t 
ISchlegel et aLllMl . 

Although the BAO signal does not evolve in linear the- 
ory, nonlinearity causes a slight shift of i ts scale and sig- 
nifica n t damping of its amplit u de (e .g. . ISeo et al.l [20081 
201(1 iSherwin fc Zaldarriagal [20121: iMcCullagh eTafl 
20131) . These effects must be modeled or corrected 
in BAO analyses. Given the unprecedented statisti- 
cal power of future surveys, it is necessary to cali- 
brate the nonlinear BAO signal accurately. Efforts 
have been made for dark matter with both TV-body 
simul ations ([Smith et al.l 120081 : IPadmanabhan fc White! 
20091) and perturbation theo ry dJeong fc Komatsul l2006t 
Crocce fc Scoccimarro: 2008; Matsubarall2008l ) . 



It is reasonable to assume that the halo bias is con- 
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stant on large scales according t o the "peak-background 
split" or "local bias" model(e.g..lFrv fc Gaztanagalll993l : 
IMo fc Whitdll^MIScheFrer fc Weinberdll998D . For pre- 
cision cosmology, however, one should quantify the scales 
of validity as well as any departure from constancy for 
halos of different masses. Investigations of the halo BAO 
signal have given hints for scale- dependence in the halo 
bias even at the BAO scale (e.g.. ISeo fc Eisensteinll2005l : 
iNoh et all 120091: iDesiacaues et all 120101) . In this work, 
we determine the halo correlation functions and power 
spectra with high precision and clearly demonstrate the 
mass-dependent halo BAO signal and behavior of the 
halo bias beyond the BAO scale. 

2. SIMULATIONS 

We run two sets of ./V-body simulations to study the 
halo BAO signal and halo bias at redshift z = 0. The ini- 
tial matter power spectrum of the first set is c alculated 
using CAMB ([Lewis, Challinor, fc Las enbv 2000) with the 
parameters: fi A =0.73, fl M = 0.27, h = 0.71, a s = 0.8, 
n s = 1, and ilb — 0.045. That o f the second set is calcu- 
lated using the fitting formula in lEisenstein fc Hul (|1998[ ) 
with the same parameters but no BAO wiggles. For con- 
venience, we refer to them as the BAO simulations and 
the no-BAO simulations, respectively. 

The simulation box is l/i _1 Gpc on each side. With 
N p = 640 3 particles, the mass resolution is 2.86 x 
10 11 h~ 1 M Q , and the softening length is se t to 50 /i~ 1 kpc . 
The simulations are run with GADGET- 2 ([Springell 12005) 
from initial redshift Zi = 100 to z = 0. 

There are 150 realizations in each set of the simula- 
tions. The initial co nditions are generated using 2LPTic 
([Crocce et al.l 120061 ) with the same random seed for a 
pair of simulations, one in each set. In this way, the ef- 
fect of sample variance can be reduced in the comparison 
between the results from the two sets of simulations. 

Halos are identified using the Friends-of-Friends group 
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Figure 1. Mass-dcpcndcnt halo BAO signal at z = 0. The sym- 
bols represent the correlation functions of M > 10 14 h~ 'Mg halos 
(filled circles), M > 10 13 /! _1 Mg halos (open circles), and dark 
matter (crosses) in the BAO simulations. The error bars show 1-tr 
errors of the mean correlation functions of 150 realizations. The 
dashed line is the linear dark matter correlation function, and the 
solid line is the nonlinear dark matter correlation function from the 
RPT fitting formula. These correlation functions are normalized 
in the range of 160-180 ft," 1 Mpc. 

finder from University of Washington's iV-Body ShopS- 
The linking length is set to 0.2 times the mean parti- 
cle separation. Given the mass resolution, we limit our 
study to halos more massive than 10 13 /i _1 M Q , which 
corresponds to about 35 simulation particles. The ob- 
tained comoving halo number density at z = is 
4.8 x l(r 4 /i 3 Mpc _3 for M > 10 13 h^M® and 2.5 x 
lCT 5 /i 3 Mpc~ 3 for M > 1O U ^ _1 M . 

3. HALO BAO AND HALO BIAS 
BAO in Halo Correlation Functions 



3.1. 



We calculate the halo two-point corr elation functions 
£ hh (r ) with the Landy-Szalay estimator (jLandv fc Szalayl 
1993) 
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where iV r is the size of a uniform random sample, Ah is 
the size of the halo sample, HH is the number of halo 
pairs with a separation between r — Ar/2 and r + Ar/2, 
HR is the number of pairs between the halo sample and 
the uniform random sample, and RR is that within the 
uniform random sample. Since the halo sample from the 
simulations has periodic boundaries and a uniform se- 
lection function, Equation (JTJ reduces, for an infinitely 
dense random sample, to 



e hh (r) = 



HH(r) V ho 



V hin (r) 



1. 



where Vbox is the volume of the simulation box, and 
Vbin(r) is the volume of the spherical shell between 
r - Ar/2 and r + Ar/2. 

The correlation functions (£ hh ) averaged over the en- 
semble of the 150 BAO simulations are shown in Fig- 
ure Q] for halos of mass M > 10 13 /i _1 M (open cir- 
cles) and M > 10 14 /i~ 1 M Q (filled circles). For com- 

1 http:/ /www-hpcc. astro, washington.edu/ 



parison, we also plot the nonlinear dark matter correla- 
tion function from the simulations, i.e., the correlation 
of all simulation particles ((£ G lm ), crosses), the linear 
dark matter correlation function (j;™ 1 ", dashed line), and 
the nonlinear dark matter correlation function from a 
convolution of £™ m with the Gaussian damping factor 
predicted by Renormalized Pertu rbation Theory (RPT, 
ICrocce fc Scoccimarrol2006l[2008l) (£™ T , solid line). Be- 
cause it is time consuming to estimate the dark matter 
correlation function from 640 3 particles with the pair- 
counting method, we assign the particles onto a density 
grid (hence the subscript G) of 200 3 cells with the cloud- 
in-cell assignment function (Hockncy & East woodlll981l) 
and then calculate the correlation function from pairs of 
the cells. 

Studies have show n that nonlinear evolution damps the 
BAO feature (e.g.. | Eisenstein et al.1 120071 : iSmith et all 
120081 : iSeo et all 120081 ). Such an effect on dark matter 
can be fairly accurately modeled with the RPT formulism 
(see, e.g., Figure[T]). Comparing the correlation functions 
from the simulations, we find that the damping effect of 
the BAO signal depends on the tracer's mass, which can- 
not be sca led away by a co n stant clustering bias. As ex- 
plained in Eis enstein et all ()2007f ). the BAO damping is 
due to the motions of matter or tracers relative to their 
initial separation, i.e., the Lagrangian displacements be- 
tween pairs. The mass-dependent damping of the halo 
BAO signal suggests that the characteristics of halo mo- 
tions are also mass-dependent. For a test, we measure the 
pairwise velocity dispersion o~ pv , which is proportional to 
the rms Lagrangian displacement in the Zel'dovich ap- 
proximation, and find that a pv = 650.1 ± 6.0kms~ 1 for 
M > 10 13 /i- 1 M Q halos and a pv = 632.5 ± 8.4km s" 1 
for M > 10 14 /i _1 M Q halos over the separations of 90- 
120 ft,~ 1 Mpc. These results appear to be consistent with 
Figurc[T]in the sense that larger pairwise random motions 
would damp the BAO signal more. However, we note 
that nonlinear evolution weakens the correspondence be- 
tweeen the pairwise velocity and the Lagrangian displace- 
ment. A more thorough investigation is needed to fully 
understand the mass-dependent BAO damping. 

3.2. Halo Bias from Correlation Functions 

The difference between the halo correlation functions 
in Figure Q] can be accommodated with scale-dependent 
biases. We estimate the halo bias in two ways: 

& CC (r)^ h V)/£ hm (r)), 
Mr)^(r))V«G m (r)>*, 



(3) 
(4) 

where £ hm is the cross correlation between the halos and 
dark matter in each simulation, and is the correlation 
function calculated from the halo density grid (analogous 
to Cg™)- The cross correlation is given by 



HM(r) V hw 
N h N p Vbin(r) 



(5) 



where HM is the number of halo-simulation particle 
pairs with a separation between r — Ar/2 and r + Ar/2. 
One might attempt to take an ensemble average of the 
ratio £g /Cg™ as an estimator for b 2 . Such a ratio can 
fluctuate wildly and even be negative in some realizations 
when the correlations are close to zero. Equation (|4| is 
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Figure 2. Halo bias from the correlation functions. Open circles 
and filled circles represent, respectively, bcc and bQ (see Equa- 
tions [3] and [4]) in the BAO simulations. Horizontal lines mark the 
average values of bQ between 50 h~ Mpc and 70 h~ Mpc. Boxes 
represent bQ in the no-BAO simulations. 

better behaved, though the covariance between £g h and 
needs to be accounted for when determining the un- 
certainties of bQ. 

Figure [2] shows the behavior of the halo bias for four 
mass ranges as labeled. Since the halo mass function 
is fairly steep at M > 10 13 /i _1 M Q , each mass range is 
dominated by halos at its low-mass end. The halo bias 
cannot be determined accurately when the correlation 
functions are close to zero. Therefore, we leave a gap 
where the errors become too large to be informative. 

For the BAO simulations, the results of the two dif- 
ferent estimators, bcc (open circles) and 6q (filled cir- 
cles), are consistent with each other over the scales and 
masses of interest. They both show a systematic off- 
set between the average halo bias over 50-70 h~ 1 Mpc 
(dashed lines) and that over the BAO peak (roughly 95- 
115 ft, -1 Mpc) for all the four mass ranges. The decre- 
ment of the bias below the transition scale depends on 
the halo mass: 4% for M > lO 14 /i -1 M and roughly 
2% for M > 10 13 /i _1 M Q . This finding is relevant to 
BAO measurements from luminous red galaxies, which 
are mostly central galaxies in halos of masses from a f ew 
times IO^/i^Mq to 1O 14 /i- 1 M (IZheng et all 120091) . 

To see whether it is a coincidence that the bias drops 
near the trough of the BAO signal, we contrast the re- 
sults from the BAO simulations (circles) with those from 
the no-BAO simulations (boxes) in Figure [5] Since the 
latter changes gradually around 90 ft _1 Mpc, we conclude 
that the drop of the halo bias in the BAO simulations is 
associated with the presence of BAO. 

The scale dependence of the halo bias around the 
BAO scale has been studied previously. For exam- 
ple, one may find a hint for a bump of the halo 
bias around 95 ft" 1 Mpc for halos with mass M > 
10 14 h^M^ in IManera. Sheth . fe Scoccimarrol (pOlof ) 
and IManera fe Gaztanagal (|2011l ). A peak of the halo 
bias (M > 1.2 x 10 14 h^Mp) nea r lOO/i^Mpc is also 
detected in lDesiacques et alT(j2010[ ). 

This work improves the statistics of the halo bias 
around the BAO scale. The mass-dependent drop of the 
bias below 90 ft -1 Mpc is well detected for halo masses 
down to 1O 13 /i _1 M . Our results are consistent with a 
flat bias over 95-115 h~ 1 Mpc. We also extend the scale 
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Figure 3. Halo bias from the power spectra. The width of each 
curve represents 1-<t uncertainties of the mean bias of the ensemble. 

of interest up to 180/i _1 Mpc where halos are anticorre- 
lated (so is dark matter). Because the correlation func- 
tions are close to zero in the range of 140-180 /i _1 Mpc, 
the uncertainties of the halo bias in this range are con- 
siderably larger than those below 120/i _1 Mpc. Never- 
theless, the results are consistent with a flat bias over 
140-180 ft" ^pc. 

3.3. Halo Bias from Power Spectra 

One can also estimate the halo bias from the halo 
power spectrum P hh (fc) and the dark matter power spec- 
trum P mm (k) with 



6 PS (fc) = ( [P hh (k)/P mm (k 



(6) 



To obtain the power spectra, we assign the ha- 
los and simulation particles to density grids of 512 3 
cells using the Daubechies D12 wavelet scaling func- 
tion, which has sup erior performan ce for measuring 
the power spectrum (jCui et al.1 120081 ). We then take 

- KoxiVh" 1 and P mm (fc) = 
where <S h (k) and <5 m (k) are 
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the discrete Fourier modes of the halo grid and the dark 
matter grid, respectively, and (. . .)g denotes an average 
within each grid. 

Figure |3] presents the halo bias determined from the 
power spectra for the same mass ranges as in Figure [H 
Although one does not expect the halo bias to be ab- 
solutely constant on large scales, it is still surprising 
that the bias of M > 10 14 /i _1 M halos falls by 11% 
from k = 0.012/i Mpc" 1 to k — 0.2 ft. Mpc -1 . For halos 
of mass M > 10 13 /i -1 Mq, the bias falls by less than 
4% over the same scales, co nsistent with the result of 
M > 1.11 x 1O 13 /i- 1 M in iPollack. Smith, fc Porcianil 
(|2012T) . 

One may notice that there are slight oscillations in the 
halo bias, which are more pronounced for more massive 
halos. This behavior is discussed in the next subsection. 

3.4. Damping of BAO and Bias Modulation 

We use the no-BAO simulations as references to exam- 
ine the damping of the BAO signal in Fourier space. The 
halo and dark matter power spectra from each simula- 
tion with BAO are divided, respectively, by those from 
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Figure 4. Upper panel: Ratio of the power spectrum from the 
BAO simulations to that from the no-BAO simulations. Filled 
circles represent dark matter, and open circles represent M > 
10 14 h — 1 Mq halos. Lower panel: Ratio of the halo bias in the BAO 
simulations to that in no- BAO simulations for M > 10 14 h~ 1 M Q 
halos. 

the corresponding no-BAO simulation. Because the pair 
of simulations share the same random seed for their ini- 
tial conditions, the sample variance is greatly reduced. 
The power spectrum ratios are shown in the upper panel 
of Figure H for M > 10 14 /i^M© halos (P hh /P$&, open 
circles) and dark matter (P mm /P^ 1 , filled circles). One 
can see that the BAO signal in the massive halo power 
spectrum suffers less damping than that in the dark mat- 
ter power spectrum. 

There is a subtle difference between Figure [Hand Fig- 
ure |U Figure [T] illustrates a mismatch between the shape 
of the halo correlation functions and that of dark mat- 
ter around the BAO scale. It does not associate the very 
presence of BAO with the mismatch. The latter could be 
an intrinsic difference between the halos and dark mat- 
ter that happens to occur around the BAO scale. Fig- 
ure [3] establishes such a link, but the uncertainty is a 
little high for M > 10 14 /i _1 M Q halos. The power spec- 
trum ratios in Figure U suppress the intrinsic difference 
between halos and dark matter and reduce uncertainties 
arising from the sample variance. It demonstrates for 
M > 1O 14 /i _1 M halos that the difference between the 
halo and dark matter BAO signals is truly due to the 
difference in the BAO damping mechanism. 

By definition, the above effect is attributed to the halo 
bias. The lower panel of Figure [4] shows the ratio of 
the halo bias in the BAO simulations to that in no-BAO 
simulations for M > 10 14 /i _1 M© halos (same as the halo 
data divided by the dark matter data in the upper panel). 
The result oscillates in phase with BAO with an ampli- 
tude of 0.5%, suggesting a modulation of the halo bias 
by BAO. This would not be sensible if halos distributed 
exactly in proportion to dark matter. However, if there 
is nonlinearity in the halo bias, as already evident from 
its scale dependence in Figures [2] and G2 the bias will de- 
pend on the shape of the dark matter power spectrum to 
some extent. 



considerably higher than that of M > 10 13 /i _1 M halos, 
which in turn is higher than that of dark matter. Since 
galaxies reside in halos, the BAO signal in the galaxy 
distribution must also depend on the mass of the host 
halos. 

We further show that the halo bias is scale-dependent 
even on very large scales. This could be a significant 
source of systematic errors when one ext racts cosmo- 
logica l information (e.g., neutrino masses, lLahav et al.l 
l2010T i . from the overall shape of the galaxy power spec- 
trum. 

Recently, iHong et all (|2012| ) detect a higher-than- 
expected BAO peak in the correlation function of galaxy 
clusters (M > 10 14 h~ 1 M Q ). Since they fix the baryon 
fraction while fitting other parameters, the excess BAO 
signal has to be compensated by a very low matter frac- 
tion. Although the difference between the halo and dark 
matter BAO signals in our study is not suffici e nt to ex- 
plain all the excess BAO signal in IHong et al.l (|2012f) . it 
should be nonetheless a significant contributing factor. 

One might expect the trend with halo mass in Figures[2] 
and [3] to continue below 10 13 h ^M^. Investigation with 
M > 1.37 x 1O 12 /i -1 M halos (|Sato fc Matsubardl20rl 
shows that the bias of such halos at z > 0.5 varies sig- 
nificantly with scales. A simple extrapolation of their 
results to z — suggests that the fractional variation of 
the bias of M > 1.37 x 10 12 /i _1 M halos could be even 
larger than that of M > 10 13 ft _1 M halos. Therefore, it 
is prudent not to extrapolate our results to lower masses. 

Besides damping the BAO signal, nonlinear evolution 
also shifts the acoustic scale sligh tly. We es t imate this 
effect using the fitting formulae in lSeo et al.l (|2008f) and 
find that the fractional shift of the BAO feature toward 
higher wavenumber is (0.28 ± 0.16)%, (0.47 ± 0.20)%, 
(0.54 ± 0.28)%, and (0.62 ± 0.23)% for dark matter, 
M > lO^h^Mq halos, M > 2 x 10 13 H^Mq halos, 
and M > 5 x 10 13 /i~ 1 M halos, r espectively . The se re- 
sults are consistent with thos e in iSeo et all (j2010T ) and 
iSherwin fc Zaldarriagal (|2012f ). The simple fitting does 
not converge well for the highest mass range in our study. 
In this case, one can use the more sophisticated fitting 
formu la for halos proposed in iPadmanabhan fc White 

Future galaxy surveys could probe effective volumes of 
a few hundred Gpc 3 , similar to the total volume of the 
simulations in this study. The effect of mass-dependent 
halo BAO signal and the scale-dependent halo bias on 
large scales could be easily detected in the galaxy statis- 
tics. In other words, it will be necessary to accurately 
calibrate these effects for precision cosmology in the fu- 
ture 
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4. DISCUSSION 

We have demonstrated with high significance that the 
BAO signal in the halo distribution is dependent on the 
halo mass. In the correlation function, the relative am- 
plitude of the BAO peak of M > 10 14 /i- 1 M Q halos is 
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